o 

o 

(N 

< 
in 



C/3 



C^ 



o 



> 

o 
o 



o 



X 



Excitation energy transfer between closely spaced multichromophoric systems: 
Effects of band mixing and intraband relaxation 
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We theoretically analyze the excitation energy transfer between two closely spaced linear molecular 
J-aggregates, whose excited states are Frenkel excitons. The aggregate with the higher (lower) 
exciton band edge energy is considered as the donor (acceptor). The celebrated theory of Forster 
resonance energy transfer (FRET), which relates the transfer rate to the overlap integral of optical 
spectra, fails in this situation. We point out that in addition to the well-known fact that the point- 
dipole approximation breaks down (enabling energy transfer between optically forbidden states), 
also the perturbative treatment of the electronic interactions between donor and acceptor system, 
which underlies the Forster approach, in general loses its validity due to overlap of the exciton 
bands. We therefore propose a nonperturbative method, in which donor and acceptor bands are 
mixed and the energy transfer is described in terms of a phonon-assisted energy relaxation process 
between the two new (renormalized) bands. The validity of the conventional perturbative approach 
is investigated by comparing to the nonperturbative one; in general this validity improves for lower 
temperature and larger distances (weaker interactions) between the aggregates. We also demonstrate 
that the interference between intraband relaxation and energy transfer renders the proper definition 
of the transfer rate and its evaluation from experiment a complicated issue, which involves the initial 
excitation condition. 

PACS numbers: 71.35.Aa; 73.63.-b; 78.67.-n 81.16.Fg 



I. INTRODUCTION 

The theory of Forster resonance energy transfer 
(FRET) between two chromophores (molecules, ions) 
with dipole-allowed optical transitions,^ and its gen- 
eralization by Dexter~ to forbidden optical transitions 
and exchange interactions between the chromophores, al- 
ready have a history of more than 50 years. This cele- 
brated theory gives an excellent description of transfer 
rates for distant chromophores with rather broad spectral 
lines. Mi^ It describes these rates in terms of the overlap 
integral of experimentally measured optical absorption 
and luminescence spectra, which makes it of great utility. 
With minor reformulations, the concept of FRET may 
also be applied successfully to the description of nonra- 
diative transitions in ions and molecules in condensed 
phases,^ as well as to energy transfer in the presence 
of a nonstationary bath relaxation^ Finally, it has been 
shown that the Forster theory also explains the efficient 
long-range energy transfer in assemblies of closely packed 
CdSe quantum dot a^i^i^" and CdSe nanocrystals assem- 
bled with molecular wiresjii systems of possible use for 
quantum computationiiSii^ 

In spite of its great success, it has been recognized since 
the 1980's that in certain situations standard FRET the- 
ory is not applicable. In particular, this holds for chro- 
mophores with narrow spectral lines and a small spectral 
overlap, such as rare-earth ions embedded in a crystalline 
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or glassy hostfi^'i^'iS^ More recently another impor- 
tant situation in which FRET theory may break down 
has been emphasized, namely energy transfer between 
two systems that both contain many interacting chro- 
mophores. This problem has drawn particular attention 
in the context of excitation energy transfer from the B800 
to the B850 ring of the photosynthetic antenna system 
];^H2.i«-^9-20-2i-22.23.24 rpj^^ g^^^ complication when deal- 
ing with systems of strongly interacting chromophores 
is that their excited states are excitons, consisting of 
a coherent superposition of the excited states of many 
molecules. Due to their spatial extent, which may easily 
exceed the separation between the two systems, the ef- 
fective interaction between an exciton state on the donor 
system and one on the acceptor system cannot be mod- 
elled as the interaction between the transition dipoles of 
both states. As a result, excitation energy transfer may 
occur from or towards a dipole forbidden (optically dark) 
exciton state, implying that a description of the energy 
transfer in terms of overlap integrals of optical spectra no 
longer holds. This breakdown of the point-dipole approx- 
imation was first pointed out by Sumi and coworkers. ^^'IS 
Treating the electronic coupling between both rings in 
LH2 as a perturbation, they derived a transfer rate be- 
tween the rings that strongly differed in magnitude from 
the Forster result and which was in good agreement with 
experiment pSSiSk 

It should be noted that LII2 is a rather special case 
of transfer between two aggregates. The reason is that 
the molecules in the B800 ring are weakly coupled to 
each other. Thus, the B800 excitations are almost 
E-mail: monomerio^i and they occur in a narrow band just above 
the upper edge of the B850 band and far (965 cm~^) away 



from the optically dominant bottom of that band>i2i2£ 
Moreover, the intermolecular interactions between both 
rings are weak, in the order of 20 cm^^iSS Given these 
special circumstances, it is not surprising that the per- 
turbative treatment of the inter-aggregate (inter- ring) in- 
teractions gives good results. 

In many cases of closely separated aggregates, however, 
a perturbative treatment of the interaction that causes 
energy transfer between them will not be valid. An in- 
teresting example is the case of nanotubular carbocya- 
nine J-aggregates that have recently been developed by 
Dahne and coworkersSSiSSiSi and which have been sug- 
gested as building blocks for synthetic light-harvesting 
systems. These aggregates consist of two walls, which 
are only a few nanometers apart. Each wall is respon- 
sible for the formation of an exciton band; even though 
their optically dominant lower edges are separated by a 
few hundred cm^^, the two bands overlap over a large 
energy range (~ 2000 cm~^)ii22. Using fluorescence and 
pump-probe experiments, fast excitation energy transfer 
between both walls has been observed. ■^^ The strong over- 
lap of the exciton bands makes it doubtful that a pertur- 
bative treatment of the interwall interactions holds for 
this example: (dark) states inside both bands may be 
close to degenerate, thus falling outside the perturba- 
tive regime. This observation also holds for the example 
of energy transfer between two linear pseudo-isocyanine 
aggregates studied by Kobayashi and coworkers i^ Like 
Sumi and coworkers, these authors focused on a break- 
down of the dipole approximation; they did not consider 
the limitations of a perturbative treatment. 

The aim of this paper is to study theoretically the en- 
ergy transfer between two molecular J-aggregates carry- 
ing Frenkel excitons. We will be inspired by the example 
of the double- wall cyanine tubes, where the bottoms of 
both exciton bands occur at different energies, but their 
central parts overlap. To keep the problem computa- 
tionally tractable, we will consider two interacting lin- 
ear J-aggregates with different bandwidths, leading to a 
crossing of both bands in their center. The chain with 
the higher (lower) band bottom is considered the donor 
(acceptor). Using this generic model, we will investigate 
the breakdown of the perturbative approach by compar- 
ing to an exact treatment, in which both aggregates form 
one exciton system and the energy transfer is associated 
with phonon-assisted relaxation within this system. We 
thus find that the crossover between the weak coupling 
(perturbative) and strong coupling (nonperturbative) sit- 
uations is determined by the separation between both 
aggregates as well as by the temperature. 

Another important issue which we address is that 
fast intraband relaxation (thermalization) between vis- 
ible and dark exciton states may obscure the observation 
of the actual energy transfer process and thereby plays a 
crucial role in the proper definition of the transfer rate as 
extracted from experiment. As a consequence, also the 
initial excitation conditions strongly affect the possibility 
to measure the transfer rate. This key role of intraband 




d 



N 



FIG. 1: The two linear J-aggregates (D and A) considered in 
this paper. Both chains contain TV chromophores, labeled n, 
which interact with each other through their transition dipoles 
(indicated by arrows). The scale of the interactions is set by 
the quantities Jd, J a, and Jda, which give the interactions 
between the various pairs of nearest neighbors. The lattice 
spacing within each chain is denoted h, while the interchain 
distance is given by d. 



relaxation in the process of energy transfer between ex- 
citonic systems seems to have gone unnoticed thus far. 



The outline of this paper is as follows. In Sec. Hll we 
first present our model of two homogeneous linear J- 
aggregates, interacting with each other as well as with 
a vibrational bath. We then introduce the general issue 
of perturbative versus nonperturbative treatment of in- 
terchain interactions by considering the various dynamic 
processes that can occur in and between weakly coupled 
or strongly mixed exciton bands. The perturbative ap- 
proach is worked out in more detail in Sec. lIIII where we 
derive the general expression for the one-phonon assisted 
transfer rate between any two exciton states located on 
different chains. In Sec. IIVI we develop the approach for 
the case of strong interchain coupling (nonperturbative 
case). Results of numerical simulations of the fluores- 
cence kinetics in both the perturbative and the nonper- 
turbative approach are presented in Sec.lVland are used 
to study the validity of the former as well as the best 
way to extract the transfer rate from experiment. Both 
the temperature and distance dependence of the transfer 
dynamics are addressed and distinction is made between 
resonant and off-resonant initial excitation of the donor. 
A comparison to standard FRET theory is made as well. 
Finally, in Sec. I VII we summarize. 



II. MODEL AND GENERAL STRATEGY 
A. Model 

We consider two parallel linear chains each consisting 
of A'^ equidistant two-level chromophores, with their tran- 
sition dipoles aligned to the chains (see Fig.QJ. The lat- 
tice spacing within each chain is denoted /i, while the 
distance between the chains is d. One of the chains will 
be referred to as the donor (D), the other as the acceptor 
(A) (see below). We will assume that both chains are 
homogeneous and we will impose periodic boundary con- 
ditions in the chain direction (the formalism may easily 
be extended to account for disorder and open boundary 
conditions). The chromophores building up the D chain 
are different from those of the A chain; in particular, we 
will assume that the transition energies and dipoles of 
the individual chromophores within the D chain all have 
the values en and /zu, respectively, while in the A chain 
they take the values ea and ha- We will account for 
the dipole-dipole interactions between all chromophores 
in both chains and also include in the model a coupling 
of the electronic excitations to a bath of vibrations. The 
latter coupling is derived from the first-order change of 
the chromophores' transition energies caused by nuclear 
displacements in the environment. In the site representa- 
tion, the resulting Hamiltonian of system and bath reads 

H = Hd + Ha + Hda 



Hhath + Hd- 



Ha- 



bath -I- ^A-bath 
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Here, Hd and Ha denote the electronic (exciton) 
Hamiltonians for the isolated donor and acceptor chains, 
respectively, while Hua is the electronic interaction be- 
tween both chains, composed of all dipole-dipole interac- 
tions between molecules of one chain and the other. In 
these terms, \n, D) (|n. A)) denotes the state in which the 
nth (n = 1, . . . , N) chromophore of the donor (acceptor) 
chain is excited and all the other chromophores are in the 
ground state. Furthermore, the hopping integrals J^, 
Jnrm and J^j^ are the various dipole-dipole interaction 
matrix elements between chromophores of the same or 
different chains. For the geometry considered here, we 
have J^,„ = -Jzj/ln-TOp, J^^ = -j^l\n-m\^, J^„ - 
Jt = 0, and J^^ = Jda [l - 2(n - mfh^/dP] / [l + (r 
m)2/i2/d2j5/2^ ^j^j^ j^ ^ 2/i|5//r\ J A = 2My/l^ and 

Jda = ^JiD^J'A/d?'■ 

Hhath describes the vibrational modes of the host, la- 
beled q and with the energy spectrum Uq (we set h^ 1). 
The operator a annihilates a vibrational quantum in 
mode q. Finally, if£)_bath and HA-hath represent the 
operators of the exciton-bath coupling of the donor and 
acceptor, respectively, where the quantities V^ and V^ 
indicate their strength. We do not provide explicit ex- 
pressions for these quantities; rather we will consider 
them on a phenomenological basis. Specifically, realiz- 
ing that in most experimental studies of aggregates the 
host is strongly disordered, we will treat these strengths 
as stochastic variables with correlation properties: 
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These relations imply that the surroundings of different 
chromophores are not correlated. 
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B. Strategy 

We now turn to a general discussion of perturbative 
versus nonperturbative approach to describe the excita- 
tion energy transfer between both chains. Throughout 
this paper, we will assume that the exciton-phonon cou- 
pling is weak compared to the intrachain dipole-dipole 
interactions, in the sense that the coherence length of the 
excitons within each chain is not limited by the coupling 
to the bath. In that case, the Bloch eigenstates of Ho 
and Ha are a good starting point for our considerations 
and the main role of the bath is to make up for energy 
differences in possible intraband energy relaxation. The 
Bloch states read {X = D,A) 
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Here, k is the wavenumber of the state, which can take 
the values 0, 1, . . . , iV — 1. 

Figure HJa) shows the donor and the acceptor bands, 
Ej^ and E^, respectively, for chains of A^ — 150 
molecules, with Ed ~ Ea (wc take this as zero of en- 
ergy) and J A — l-14Ji3. Only half of the Brillouin 
zone is plotted, as the dispersion relation is symmetric 
around k = 7V/2. We see two exciton bands that cross at 
their center. The bottom of each band occurs at A: = 0, 
which is the superradiant state that contains all oscilla- 
tor strength to the ground state; all other exciton states 
are optically forbidden. Thus, as long as the interactions 
between the chains are weak, the fluorescence spectrum 
consists if two J-bands, whose positions are indicated by 
the dots on the fc = axis. By definition, we choose our 
labeling such that E^^^ > E^^^. The energy separation 
between both peaks will be denoted A and is (for N ^ 1) 
given hy A — Ed — £a + 2({3){Ja — Jd), with ({x) the 
Riemann zeta function. 

The two chains are coupled electronically by the inter- 
action HuA (Eq- if^ ). which on the basis of Bloch states 
is diagonal, i.e., Hda — J2k Jkk^\^'D){^i^\ +h.c., with 
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In general, this expression must be evaluated numerically 
(see Fig. ^c) for an example) . Only in the limiting (and 
rather unphysical) case of d <C /i it is easily seen that 
Jkh — '^DA = IJ-DiJ'A/d? , independent of k. 

As long as the interactions J^^ are weak (large sepa- 
ration d), it seems reasonable to apply the usual pertur- 
bative approach to them, using the bands of the isolated 
donor and acceptor chains as starting point. We will fol- 
low this route in Sec. IIIII where we also will treat the 
exciton-phonon interactions perturbatively. In this ap- 
proach, the rate of energy transfer between an arbitrary 
donor state and acceptor state is calculated using second 
order perturbation theory, involving two steps: (i) The 
excitation is transferred from the donor to the acceptor, 
(ii) A phonon-induced scattering occurs within the accep- 
tor band (scattering within the donor followed by transfer 
to the acceptor is possible as well). The overall process 
should conserve the total energy. Both steps are schemat- 
ically indicated in Fig. [5Jb), in which we zoomed in on 
the small- A; part of the bands given in Fig.j^l^a). As in our 
example the interchain interaction is diagonal in fc, the 
transfer step is a vertical transition between both bands. 
We notice that aside from energy transfer between both 
chains, the interactions with the phonon bath also give 
rise to relaxation of the exciton states within both bands; 
details of this process are studied in Ref. ISfl. 



If we focus on the optically dominant bottom states 
of both bands, it seems that the perturbative treat- 
ment of the interchain interactions is valid as long as 
l'^fc^fc=ol "^ ^- Generally speaking, this criterion is 
not sufficient, however. Depending on temperature and 
initial excitation condition, the transfer process may in- 
volve the optically forbidden states higher in the exciton 
bands. As near the band center the donor and accep- 
tor states get arbitrarily close in energy, a perturbative 
approach necessarily fails there. For this situation, one 
has to resort to a nonperturbative treatment, in which 
both bands are mixed. As long as \JJ?^q ^.^qI ^ ^' ^^^ 
amount of mixing will be weak for the superradiant band 
bottoms, but strong near the center. 

For our case of translational symmetry, the mixed 
eigenstates of Hd + EIa + EIda read (see, e.g., Ref. i37f) 
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with energies 

Ei = \{E^ + Et)±\{E^^Et).JTTVk- (7b) 

Here, ± label the new (decoupled) exciton bands and 
■q^ = 4| J^^|2/(^f - E^f is the quantity that charac- 
terizes the amount of mixing at wavenumber fc. The new 
bands are plotted in Fig.j^fc) for the same parameters as 
in Fig. |2Ia), accounting for an interchain interaction of 
strength JoA = 0.535 J^) (which corresponds to a small 
interchain separation of d = h). The coupling J^ is 
depicted as well. As we see, at the band bottoms, the 
coupling is very small compared to A and the states are 
hardly mixed {rik=a = 0.2). Hence, we may still refer to 
the fluorescence coming from the bottom of the upper 
(lower) band as the donor (acceptor) fluorescence. At 
the center of the bands, the mixing is very strong and 
a band anticrossing occurs. Also after mixing, the only 
states with oscillator strength occur at the band bottoms. 

After accounting for the band mixing, the only remain- 
ing dynamics that can take place is phonon-assisted scat- 
tering of the new exciton states, leading to intraband and 
interband relaxation of the excitation energy (see pro- 
cesses indicated in Fig. |3d)). In Sec. II VI this approach 
will be further specified, using the Fermi golden rule to 
account for the exciton-phonon interaction. 

To end this section, we stress that the special condi- 
tions imposed in our model, such as chains of equal length 
and periodicity, periodic boundary conditions, and the 
absence of disorder in the electronic part of the Hamil- 
tonian, may be relaxed without affecting the above for- 
malism and the distinction between the perturbative and 
nonperturbative approach. In fact, in all expressions pre- 
sented in Sees, mil and llVI we will use a general notation 
for the exciton wave functions, V'foiJ Vkn' ^^'^ (f>vn, so 
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FIG. 2: (a) Exciton bands for donor (solid) and acceptor (dashed) chains oi N — 150 molecules, in the absence of interchain 
coupling for the case Ja = 1.14Jd and setting ea = sd = 0. Note that we have only plotted half of the Brillouin zone, as it 
is symmetric around k = N/2. The dots at fc = indicate the positions of the absorption peaks of donor and acceptor band, 
respectively, (b) Expansion of the small-fc part of panel (a), with schematic indication of the two-step perturbative view of 
the interchain energy transfer described in the text, (c) Exciton bands for the same system as in panel (a), but now after 
accounting for the mixing between the eigenstates of both chains due to the interchain interaction with strength Jda = 0.535 Jo 
(corresponding to an interchain separation d — h.) The fc-dependence of the interchain coupling Jj.j, is depicted as well, (d) 
Expansion of the small-fc part of panel (c), with schematic indication of the various phonon-assisted intraband and interband 
relaxation processes that together are responsible for the exciton dynamics within the nonperturbative picture. 



that these expressions keep their validity under general- 
ized conditions, as long as the proper exciton eigenfunc- 
tions of the generalized Hd, Ha, and Hd + Ha + Ho a 
are used as input. In that case, k and v refer to the 
appropriate quantum numbers, which do not necessarily 
have the meaning of quasi-momenta. 



and the exciton-bath interactions take the form 
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Hda^ J2 Jkk^\k,D){k',A\+h.c., (8a) 
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III. PERTURBATIVE APPROACH 



In this section, we follow the perturbative approach 
outlined in Sec. IIIBI In this approach, the Hamiltonian 
of the unperturbed system reads Hq = H d + H a + H\,a.t\i, 
while H' = Hda + ^D-bath + ^A-bath represents the 
perturbation that induces transitions between the eigen- 
states of Hq. On the basis of exciton eigenstates of the 
noninteracting chains (Eq. ^, the interchain coupling 
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E T,^kl,\k,X){k',X\a, + h.c., (8b) 
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with V^lg = E„ 14^^fe^'P^„ and X^D,A. 

We aim to consider the energy transfer between any 
donor state |fc, D) and any acceptor state \k' , A). Because 
the transfer operator H' does not produce such transi- 
tions within first-order perturbation theory, we have to 
resort to the second-order term. The corresponding ex- 



pression for the transfer rate reads 
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Here, \i) = \k,D)\{nq},) and |/) = |fc',A)|{nJ/) are 
the initial and final states, respectively, where {riq} de- 
notes the set of occupation numbers of the vibrational 
modes. Ei — E^ + fi^ and Ej — E^, + Vtf are the corre- 
sponding energies, with fij and 17/ denoting the energies 
of the bath in the initial and final states, respectively. 
Furthermore, s labels the intermediate states |s), with 
energies Eg and the quantity p{Ei) is the equilibrium 
density matrix of the bath's initial state. Finally, the an- 
gular brackets indicate that we average over the stochas- 
tic realizations of the surroundings of donor and acceptor 
chromophorcs. 

Evaluating the expressions in Eq. and accounting 
for the stochastic properties Eqs. ||2J) of the exciton-bath 
couplings, we obtain 
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where ^{lo) is the vibration spectral density of the bath, 
which we assume to be identical for the donor and accep- 
tor. It is given by 



V{uj) = 27r 
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where X = D^A. The quantity O^, denotes the prob- 
ability overlap of the donor states \k,D) and \k' ,D) for 
X = D and of the acceptor states \k,A) and \k' ,A) for 
X^A: 
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Finally, n{ujq) — [exp(ajg/r) — 1]^-'^ is the mean occupa- 
tion number of the vibrational mode q (the Boltzmann 
constant fcs = 1). 

The transfer rate Eq. (|l()(l reflects the two-step pro- 
cesses introduced in Sec. Ill Bl The first term corresponds 
to the process in which the initial exciton of wavenumber 
k on the donor chain is scattered into the exciton state 
k" on the donor under the creation or annihilation of a 
phonon q, followed by the transfer of the exciton k" to 



an exciton k' on the acceptor chain. Likewise, the second 
term in Eq. (|10|l derives from the process in which the 
donor's initial exciton k is first transferred to the accep- 
tor chain fc", followed by a phonon- induced scattering to 
the final acceptor state k'. 

Equation H1Q(I clearly indicates that in general the en- 
ergy transfer may occur between any pair of donor and 
acceptor states, independently of whether those states 
are optically allowed or forbidden: neither the prob- 
ability overlap O^f., nor the transfer interactions J^'^ 
vanish for general k and k'. Equation (|10|l reduces to 
Forster's formula only if the distance d between the 
chains is larger that their lengths Nh. Indeed, in that 
case the interchain interactions are well-approximated 

by Jg^ = ifTDf^A/d^) Eti Vl^^ ELi ^t',n. Where the 
sums correspond to the dimensionless transition dipole 
moments of the states \k, D) and \k', A) . For our ordered 
chains, these dipole moments are giant for k — k' = 0, 
while they vanish for all other states. Thus, the term 
Jk=o k'=o ^^ dominant. The energy transfer between the 
chains may then be viewed as Forster-type due to the 
overlap of the donor's side-band fluorescence and the 
zero-phonon acceptor absorption and vice versa. 

The importance of the forbidden exciton states in the 
energy transfer between excitonic systems was mentioned 
for the first time by Sumi and co-workersi^iiS and later on 
received attention from several authors iSii2iiS2iS2iS4i21 We 
note that in Refs. [l^ and [l^ the problem was treated in 
a more general way than we did above: only the transfer 
interactions J^„ were considered as perturbations, while 
the exciton-bath couplings were accounted for by a self- 
consistent second-order evaluation of the self-energy. We 
do not use such a self-consistent treatment of the exciton- 
bath interaction in this paper, as it turns out impossible 
to combine it with the nonperturbative treatment of the 
transfer interactions which we will investigate in the next 
section. 

The expression for Wk'A,kD obtained above allows one 
to calculate the energy transfer rate between any pair of 
donor and acceptor states. From the occurrence of the 
energy denominators in Eq. H10|l. however, it is clear that 
the perturbative approach of the interchain interactions 
is restricted to situations where the donor and acceptor 
bands do not cross. A band crossing or (for disordered 
systems) a band overlap will lead to small denominators, 
breaking the result of perturbation theory and giving rise 
to rates that are large compared to the energy spacing of 
the states involved. This motivates the nonperturbative 
approach presented in the next section. 



IV. NONPERTURBATIVE APPROACH 

In this section, we work out in more detail the non- 
perturbative approach outlined in Sec. IIIBI in which the 
interchain interactions are accounted for exactly, through 
mixing the exciton bands of the isolated chains. We will 
label the 2A^ mixed eigenstates oi Hex = Hu + Ha + Hda 



with a greek index v = 1, . . . , 2N; they take the form 



2N 



(13) 



where |n) = |«, -D) for n = 1, . . . , TV and \n) = \n — N, A) 
ior n = N + 1, . . . , 2N. The corresponding energy is de- 
noted E^. For the special highly symmetric model con- 
sidered in this paper, the explicit forms of the eigenstates 
and energies have been given in Sec. Ill Bl ahead v. where 
the v index should be identified with the kzt labeling. In 
the current section, the notation is kept general. 

In the new representation, the exciton-bath coupling, 
which is the perturbation that causes the exciton dynam- 
ics, reads 

2N 

HD-hath+HA-hath = ^ ^Vf,„q\fj.){i'\aq+h.c., (14) 

fi^l'—l q 



qVfj.n'riJn, 



where, as before, the cou- 



With Vf_,^q = Z]„=l ^n 

pling strength Vnq is considered a stochastic function of 
the site index n with properties similar to those given by 



Eqs. ®: (Vnq) =0 and 



'^ mq 



V* ) 

nq/ 



IKI 



As discussed in Sec. IIIBI the dynamics of the excitons 
in the new representation is caused by their scattering on 
phonons. In order to describe this process, we will use a 
Pauh Master Equation for the populations Pv{t) of the 
exciton states: 



P. = -i.P. + V {w,^p^ - w^,p,) 



M 



(15) 



Here, 7,^ — ^qF^, is the radiative decay rate of the state 
\v)^ where 70 is this rate for a single chromophore (for 
simplicity taken identical for donor and acceptor) and 



F^ = \Yl,n=i^vn) denotes the dimensionless oscillator 
strength of the exciton. Furthermore, the scattering rates 
W^u are obtained using Fermi's golden rule, taking into 
account the stochastic properties of the exciton-bath cou- 
pling strength Vqn- They are given by^S^ 

W^, = V{\E^~E,\)Of,, 

\ + n{E,-E^) , E,>E^, 
X < (16) 

n{E^~E,), E^>E,, 

where I'(w) is the spectral density of the bath, given by 
Eq. Hll|) with V^ replaced by Vq. Furthermore, O^i^ is 
the probability overlap of states ^ and i^, 

2N 



Of,^ = ^ \(f>f,n\'^\(l>u 



(17) 



n=l 



The explicit expression for O^j, in the case of homoge- 
neous chains of equal length is given in the Appendix. 



To end this section, we note that in general the best 
way to probe the excitation energy transfer is to follow 
the fluorescence kinetics of donor and acceptor. Even if 
we allow for inhomogeneity, two coupled aggregates with 
a clear separation between their individual J-bands will 
lead to a mixed system that still has two J-bands, possi- 
bly renormalized in magnitude and shifted somewhat (see 
Fig. [21 for the special case of ordered chains). In terms 
of the energies of the mixed eigenstates and the solution 
to the Pauli master equation, the time-dependent fluo- 
rescence spectrum reads 



iiE,t) = (Y,7.p.imE 



(18) 



where the brackets denote the average over disorder, if 
present. One may then define the total donor (acceptor) 
fluorescence as the integral over the highest (lowest) peak 
in the spectrum. These two quantities will be analyzed 
in Sec. [V] as a function of time following some initial 
excitation of the donor chain. 



V. RESULTS AND DISCUSSION 

We now apply the formalism developed in the previ- 
ous sections to study the energy transfer rates and flu- 
orescence kinetics for the model introduced in Sec. Ej 
namely two parallel homogeneous chains of N chro- 
mophores, with periodic boundary conditions. In all 
examples, we will choose N = 150, Ea = £d and 
J A — 1.14J/3, which implies that 11 a — l-07/i£); the 
value of Jda depends on the interchain distance d, which 
we will vary. Given the above parameters, we have 
Jda = {iiA/nD)ih/dfJD/2 = Q.hih{h/dfJD. 

For the spectral density of the bath we choose 



UJ 

D{uj) — Wo — exp 

UJc 



UJ 
UJc 



(19) 



which represents a model function with Ohmic (i.e., 
linear) behavior?S for frequencies up to a cut-off fre- 
quency UJc- The overall prefactor Wq is a measure 
of the exciton scattering amplitude imposed by the 
phonons. A spectral density similar to Eq. H19|) . with 
UJc in the order of 100 cm~^, has been used success- 
fully to fit the optical dynamics in photosynthetic an- 
tenna complexes i^242i^i*^ Also, Eq. H19|l without a cut- 
off has been used to explain the optical dynamics of ag- 
gregates of the dye 3,3'-bis(sulfopropyl)-5,5'-dichloro-9- 
ethylthiacarbo-cyanine (THIATS) measured between 
and 100 Ki^ In the remainder of this paper we will use 
UJc = 0.2 Jo, which is reasonable for J-aggregates, where 
Jd typically lies in the range 500-1000 cm~^. Wq and 
Jo will be kept arbitrary, where Wq will serve as unit 
of time and Jjj as unit of temperature. For example, for 
Jd = 800 cm"-'^, the choice T/Jd — 0.25 agrees with 
room temperature. 



The above specifies all input necessary to determine 
the transfer rate Wk'A,kD [Eq- CDIl] between arbitrary 
states in the perturbative approach, as well as the relax- 
ation rates W^i, [Eq. (|16|l ] in the nonperturbative treat- 
ment. In the latter case, the transfer from the donor to 
the acceptor involves the relaxation rate Wk'-.k+- 



A. Nonperturbative-to-perturbative crossover 

In order to assess the validity of the perturbative ap- 
proach, we calculate the total transfer rate from donor 
to acceptor chain within both approaches, assuming that 
the initial exciton populations of the donor manifold are 
in thermal equilibrium, i.e., we assume that the intraband 
relaxation is fast compared to the energy transfer process. 
For the nonperturbative approach, this means that we 
use as initial condition Pk+{t = 0) = Z^^ exp{—Ek+/T), 
where Z+ — ^j. exp(—Ek+/T), while the acceptor band 
is not populated, Pk-[t = 0) = 0. Then, the effective 
quantity 



W-+ = 



1 

zl 



2_^exp 



E, 



k+ 



T 



Wk' 



,k+ 



(20) 
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FIG. 3: Donor to acceptor energy transfer rate in the non- 
perturbative approach (thick lines) and the perturbative ap- 
proach (thin lines), as a function of the interchain distance d 

for various temperatures T. These rates are given by W ^ 

[Eq. H2()|l ] and Wad [Eq. I|21|l ]. respectively. All parameters 
as given at the beginning of Sec.lVI 



may be associated with the energy transfer rate from the 
donor to the acceptor band at the initial stage of the 
transfer, when the back-transfer is negligible. In the per- 
turbative treatment, the analog of this effective transfer 
rate is given by 



k,k' ^ 



Wk'A,kD 



(21) 



with Zd = Z]fe6xp(— ii^^/T). By construction, W-+ -^ 
Wad in the limit of large interchain distance d, i.e., small 
interchain interactions {rjk ^1). 

Figure 13 shows the results of our calculations for W ^ 

(thick lines) and Wad (thin lines) as a function of the in- 
terchain distance d. From this figure it is clearly seen that 
the perturbative approach strongly overestimates the en- 
ergy transfer rate for small d, as expected. It is also 
seen that in the limit of large d, both treatments indeed 
are in perfect agreement. However, the value of d be- 
low which the perturbative approach fails, increases for 
growing temperature. This finds its natural explanation 
in the band mixing. Let us first consider the situation at 
T = 0. Then only relaxation from the donor band edge 
state |0-f ) to states in the acceptor band below the donor 
band edge contribute to W |-. For our choice of param- 
eters, even at distances as small as d = h, these states 
in the vicinity of the band edges are only weakly coupled 
iVk^o is rather small, see Sec. IIIB|l . This explains why 
the difference between the perturbative and nonpertur- 
bative results quickly vanishes for d > h. This difference 
only is considerable for d < h, because the interchain 
interactions then quickly grow (~ 1/d^), increasing the 
band mixing even for the band edge states. 
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FIG. 4: Energy transfer rate as a function of interchain dis- 
tance at T = according to Forster's formula (dashed line, 
see text for definition) compared to the exact result [solid, 
Eq. 1201 1 and the perturbative one [dashed, Eq. 12H l. All pa- 
rameters as given at the beginning of Sec. [3 The inset shows 
the convergence to Forster's result for very large distances. 



Upon increasing the temperature, the initial popula- 
tion will spread to higher k states in the donor band. 
For these states, the band mixing becomes increasingly 
more important, even for growing values of the distance 
d. This effect finds its origin in the fact that the en- 
ergy separation between the isolated donor and acceptor 
bands decreases with growing k, while the interchain in- 
teraction J^jP grows (cf. Fig. 121). This explains why the 



range of d over which the perturbative approach fails, 
increases with growing temperature. 

It is of interest also to compare the above derived en- 
ergy transfer rates with Forster's formula, which is ob- 
tained from Eq. (|21|) under the assumption that J^"^ — 
{pLDtJ'A/d^)N5}.o (cf. Sec. IIII|I . Here, we hmit ourselves 
to the rate at T = 0, for which Forster's rate reduces to 
W^jj = 1911Wo(/i/d)6. This result is plotted as a func- 
tion of d in Fig.0I together with W-+ and Wad aXT — Q. 
Clearly, the Forster result gives an enormous overestima- 
tion of the exact as well as the perturbative rate for in- 
terchain distances smaller than roughly the chain length. 
This is not surprising, because at distances small com- 
pared to the chain size, the Jj?j^ defined below Eq. (jH3l, 
does not reduce to the interaction between the excitons' 
transition dipoles. 



B. Fluorescence kinetics 

In the above, we characterized the excitation energy 
transfer by a single rate, which is possible if we assume 
fast equilibration inside the donor and the acceptor band, 
and neglect the back transfer from the acceptor to the 
donor. In general, one cannot rely on these assumptions. 
As we mentioned in Sec. lIIBI falso see Sec. lIV|) . the most 
straightforward way to characterize the energy transfer 
between donor and acceptor in experiment is to follow the 
fluorescence kinetics of both subsystems. In this section, 
we will use this approach. In all cases, the calculations 
were done using the nonperturbative method, i.e., ac- 
counting for band mixing and for relaxation within and 
between the mixed bands. We will not assume a priori 
that equilibration takes place on a timescale shorter than 
the one for energy transfer. As only the bottom (fc = 0) 
states of the mixed bands have oscillator strength, both 
the absorption and the fluorescence spectrum consist of 
a single i5-peak. Thus, the total donor and acceptor flu- 
orescence intensity as a function of time are given by 
loit) = ^o+Po+it) and /^(i) = 7o-Po-(Oi respectively. 
These are the quantities that will be analyzed in the fol- 
lowing and related to the intraband relaxation and ex- 
catitation energy transfer. 



1. Resonant excitation 

We first study the fluorescence kinetics for resonant 
excitation, i.e., assunung that initially the superradiant 
state |0-|-) of the donor band is excited by the pump. 
Following this excitation, the population of this state has 
two channels to relax, namely scattering to higher states 
in the donor band (intraband relaxation) or scattering 
to the acceptor band (energy transfer). Once transferred 
to the acceptor band, the excitation may undergo re- 
laxation within that band, or it may transfer back to 
the donor. Figure shows the donor fluorescence kinet- 
ics /_D(i) resulting from the interplay of these processes 



at different temperatures and interchain distances, calcu- 
lated by solving Eq. ((T3|l with the above initial condition. 
In doing so, we neglected the radiative decay rates 7^ of 
the exciton states, assuming this relaxation channel to be 
much slower than all the others. As a result, at nonzero 
temperatures, Po+(i) eventually reaches a flnite value, 
which is in accordance with the Boltzmann equilibrium 
over donor and acceptor states. We notice that even the 
neglect of the superradiant emission rate for chains of the 
order of 100 molecules (10— 100 ps lifetime) does not limit 
the validity of the results presented in this section. Typ- 
ical values for I/Wq are in the 1 — 10 fs time scalej^ so 
that at all times considered in Figs.[5land|Hlthe radiative 
decay indeed has a negligible effect. 

In Fig. 13 (b), (c), and (d) one can distinguish two stages 
in the donor fluorescence kinetics, prior to reaching the 
Boltzmann equilibrium. These stages are characterized 
by different time scales. The first stage, indicated by I 
in Fig.[SJc), is a fast intraband relaxation. Here, the ini- 
tially created population of the superradiant donor state 
|0, -I-) rapidly scatters to the higher-lying dark states 
|fc, -I-) , resulting in a fast decay of the donor's fluorescence 
intensity. Obviously, this process gets more pronounced 
with increasing temperature. The second stage, indicated 
by II in Fig.|3Jc), is noticeably slower. This stage is asso- 
ciated with the interband relaxation, i.e., with the energy 
transfer between the donor and acceptor bands. We note 
that at T = the upward intraband relaxation is absent 
and the monoexponential fluorescence decay directly re- 
flects the transfer rate calculated through Eq. (|20l) at 
T = 0. 

As is seen, stages I and II are easily separated as long 
as the interchain distance d > /i; it is then possible to de- 
termine a meaningful effective energy transfer rate from 
the stage II decay of the donor fluorescence. At d = /i 
[Fig. [Sja)], it is hardly possible to distinguish the two 
stages; the fluorescence kinetics only exhibits a fast initial 
drop, almost directly followed by the Boltzmann plateau. 
The reason is that for d < h, the bare donor and acceptor 
bands are strongly mixed (see Fig. 13). As a result, the 
intra- and interband relaxation rates are of the same or- 
der, making it impossible to distinguish their signatures 
in the fluorescence kinetics. We also observe that the 
intraband relaxation is much more sensitive to changing 
the temperature than the energy transfer process. This 
is due to the fact that this relaxation is directly sensitive 
to the number of thermally accessible dark states in the 
donor band. 

The above statement that the stage II part of the donor 
kinetics indeed yields a meaningful measure for the en- 
ergy transfer rate to the acceptor, is corroborated by 
Fig. El in which for three different separations at two 
different temperatures the donor fluorescence kinetics is 
plotted again (thick solid curves), and is compared to 
the monoexponential decay obtained through Eq. (|20|l at 
the same temperature (dashed line). In Figs. |nib,c,e,f) 
stage II can be distinguished and this part of the kinetics 
is (roughly) parallel to the monoexponential curve. The 
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FIG. 5: Kinetics of the donor fluorescence loii) following resonant excitation [Pk+{Q) = 5ko] for different temperatures T and 
distances d between the chains, as indicated in the panels. All other parameters as given at the beginning of Sec. EI The thick 
line in panel (c) separates the two kinetic stages (I and II) distinguished in the text. 



agreement deteriorates with decreasing separation d and 
with growing temperature. In panel (e) stage II already 
is rather hard to see and the agreement with the mo- 
noexponential curve is only fair. In panels (a) and (d), 
stage II cannot be distinguished, and no part of the fluo- 
rescence kinetics shows a reasonable agreement with the 
monoexponential decay line. 

From the above it may be concluded that in the case 
of strong interchain coupling, under conditions of reso- 
nant excitation, the kinetics of the donor fluorescence is 
not a good tool to determine the energy transfer rate. 
One may expect that under these conditions, the accep- 
tor fluorescence yields a better tool, because this quan- 
tity is influenced less directly by the upward relaxation 
within the donor band. For this reason, we also plot- 
ted in Fig. the acceptor fluorescence as a function of 
time (thick dash-dotted curves). These curves represent 
7o-[Po-(^ -^ oo) — Po-it)], where Po^{t — + oo) equals the 
thermal equilibrium value. We observe that even in the 
case of strong coupling and high temperature, the accep- 
tor fluorescence indeed seems to be affected less by the 
initial thermal relaxation part and more directly reflects 



the energy transfer rate as calculated assuming equilibra- 
tion (dashed lines). 



2. Off-resonant excitation 

We next turn to the case of off-resonant excitation. 
To this end we also plotted in Fig. (thin sohd curve) 
the donor fluorescence after the system has been brought 
(artificially) in the initial state with Pk+{t) — Sk^N/2 ^^nd 
Pk-{t) — 0, i.e., the center state of the donor band is ex- 
cited. As for the case of resonant excitation, we observe 
that at low temperatures and (or) large interchain sep- 
arations, the kinetics exhibits two stages prior to reach- 
ing the Boltzmann plateau. The first stage, intraband 
relaxation, now manifests itself through the increase of 
the donor fluorescence from its initial value zero, because 
the donor population needs time to reach the superradi- 
ant bottom state of the donor band. We note that, in 
spite of the many relaxation steps needed to reach the 
bottom, for low temperatures this first stage occurs on 
the same timescale or even faster than for the case of res- 
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FIG. 6: Kinetics of the donor fluorescence for resonant excitation (thick solid curve) and off-resonant excitation (thin solid 
curve) calculated for various distances (d//i = 1,2 and 5 from left to right) and temperatures (T = 0.05Jd and 0.25Jd for top 
and bottom rows, respectively). All other parameters as given at the beginning of Sec. [3 The thick dash-dotted curve shows 
the difference between the acceptor's fluorescence intensity and its thermal equilibrium value (see text). Finally, the straight 
lines represent the donor decay assuming fast equilibration over the donor band in the nonperturbative case [dashed curve, 
Eq. I20II ] and the perturbative limit [dotted curve, Eq. 12H 1. 



onant excitation. As a result the stage II process, related 
to the energy transfer, reflects the same timescale as for 
resonant excitation. 

By contrast, at higher temperatures the distinction be- 
tween the first and second stage of the kinetics is less 
clear: there is a more gradual change between the two, 
which also creates the impression of a longer timescale for 
the intraband relaxation. We explain this as follows. As 
near the band center the mixing of donor and acceptor 
bands is very strong, the intraband and interband relax- 
ation there happen at the same timescale. As a result, 
off-resonant excitation is followed by very fast distribu- 
tion of the excitation over donor and acceptor bands, 
in both of which the population relaxes to the bottom. 
Thus, as opposed to the case of resonant excitation, the 
energy transfer problem at the bottom of the band takes 
place in a situation where a large part of the population 
already has ended up in the acceptor band. In this situ- 
ation, the back-transfer from the acceptor to the donor, 
which takes place at temperatures of the order of the 
band edge detuning A, feeds the donor fluorescence af- 
ter the initial relaxation, which prolongs the time during 
which the donor fluorescence increases and slows down its 
overall kinetics. Comparing to the dashed lines in Fig.|Sl 
we observe that at the higher temperatures no part of 
the kinetics curve for off-resonant excitation reflects the 



energy transfer rate calculated through Eq. H2U|I . Thus, 
at high temperatures (comparable to A or higher) off- 
resonantly excited donor fluorescence seems not to be a 
good ruler for the energy transfer rate. At low temper- 
atures, this technique is more useful and then yields the 
same transfer rate as donor or acceptor fluorescence after 
resonant excitation of the donor. 

Without showing details, we mention that in the case 
of off-resonant excitation also the kinetics of the acceptor 
fluorescence suffers from the interplay between intraband 
and interband transfer and therefore also does not yield 
a proper tool for measuring the energy transfer rate. 



VI. SUMMARY AND CONCLUSIONS 

We theoretically studied the excitation energy trans- 
fer between two linear chains of dye molecules, whose 
primary excitations are Frenkel excitons. The superra- 
diant band bottoms of both chains were detuned by a 
value A; the chain with the higher (lower) band bot- 
tom was considered the energy donor (acceptor). In ad- 
dition to electronic interchain interactions, our theory 
accounts for weak exciton-phonon coupling. Two meth- 
ods were investigated. In the perturbative method, we 
used the common approximation of treating the inter- 
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chain interactions as a perturbation, which also is the 
basis for the standard FRET formalism. In the nonper- 
turbative approach, we allowed for mixing of the donor 
and acceptor exciton bands and considered the various 
phonon-induced relaxation mechanisms in and between 
these bands to describe the energy transfer process. 

Using the perturbative approach, we have confirmed 
that energy transfer between multichromophoric sys- 
tems can occur from and towards dipole forbidden 
statesi^^i^^i^^i^^i^^i^'^i^'^i'^'^ We have focused, however, on 
two other aspects specific to multichromophoric systems, 
namely the effects of band mixing and intraband relax- 
ation on the energy transfer. The first has been studied 
by comparison of the perturbative and the nonperterba- 
tive method. Not surprisingly, we find that the perturba- 
tive approach looses its validity when the interchain dis- 
tance decreases (and the interchain interactions grow). 
The perturbative result for the energy transfer rate then 
strongly overestimates the exact one that accounts for 
band mixing. The distance at which the perturbation 
theory breaks down increases with temperature, because 
the effect of band mixing is larger for the states inside the 
bands, which get populated upon heating. We also com- 
pared our results to standard FRET theory, which turns 
out to lose its validity as soon as the interchain separa- 
tion falls below the chain length (or exciton delocalization 
length, in case one considers disordered chains). 

From our study it turns out that the process of energy 
transfer between multichromophoric systems is strongly 
affected by intraband relaxation. This is most obvious in 
the kinetics of the donor and acceptor fluorescence after 
pulsed excitation. As these quantities are natural choices 
for probing energy transfer, it is important to realize that 
this kinetics does not only reflect the actual transfer pro- 
cess of interest, but also is affected by thermal relaxation 
inside the donor and acceptor bands. Only if both pro- 
cesses can be separated, due to different time scales, one 
may extract the energy transfer rate from such exper- 
iments. Our results suggest that the best way of de- 
termining the transfer rate between two J-aggregates is 
to measure the fluorescence kinetics of the acceptor's J- 
band after resonant excitation of the donor band. Under 
these conditions, the approach of the acceptor's fluores- 
cence intensity grows towards its equilibrium value al- 
most mono-exponentially, allowing for a meaningful def- 
inition of the transfer rate. We also found that the thus 
extracted rate agrees well with the transfer rate obtained 
when assuming equilibration over the donor band. If, due 
to relaxation to the ground state, measuring the accep- 
tor's equilibrium fluorescence intensity turns out to be a 
problem, one may also resort to analyzing the decay of 
the time-derivative of this intensity. The first component 
of that decay should then refiect the transfer rate, while 
the second reflects the relaxation to the ground state. 

We note that throughout this paper we have neglected 
the homogeneous linewidths F^ of the individual exci- 
ton transitions. Within our model, these widths result 



from exciton-phonon scattering.^- It may be estimated 
that within the context of excitation energy transfer these 
widths may indeed be neglected (i.e., F^ <C J^/;"^) as long 
as Wo <^ 2ttJ'^/uJc- 

The results of this paper were derived using the sim- 
plest possible model of two interacting homogeneous 
chains with periodic boundary conditions. The same is- 
sues of band mixing and intraband relaxation will play 
an important role for more general systems, where one 
allows for energy and interaction disorder and (or) for dif- 
ferent geometries. In fact, all expressions in Sees. IIIII and 
IIVI also hold for such more general situations, provided 
that the various wave functions (/j£j, (/j^„, or (f)^n as well 
as the energies E^ , E^, or E^^ are replaced by the eigen- 
vectors and eigenenergies of the corresponding exciton 
Hamiltonian. An example of great interest to analyze 
next, is a pair of concentric cylindrical aggregates. As 
explained in the Introduction, our model of two parallel 
linear aggregates was inspired on this experimentally re- 
alized situationi2Si2i Varying the distance between both 
walls in this system is possible by altering molecular side 
groups and (or) changing the solvent. 
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APPENDIX A: SOME EXPLICIT EXPRESSIONS 
FOR HOMOGENEOUS CHAINS 

In this Appendix, we present analytical expressions for 
the probability overlap functions valid for the ordered 
case. 

Substituting the eigenfunctions Eq. Q) into Eq. I|12f) . 
the probability overlap O^, is obtained as 






1 

TV 



(Al) 



Likewise, from Eq. (|7a|l and Eq. H17|) we obtain the prob- 
ability overlap O^^, [with /i and v chosen from (fc, -I-) or 
(k, — )] in the form 



v={k',±), 



Ou. 



1 

2N 



^{^ + rik){l + r]k') 



^J{l + rik){l + r}k') 



M= (^,±) , 



z.= (fc',T) 



(A2) 



13 



* On leave from S.I. Vavilov State Optical Institute, Birzhe- 
vaya Liniya 12, 199034 Saint-Petersburg, Russia. 

^ Th. Forster, Ann. Phys. 2, 55 (1948); in Modem Quantum 
Chemistry, Pt. Ill, edited by O. Sinanoglu (Academic, New 
York, 1965). 

^ D. L. Dexter, J. Chem. Phys. 21, 836 (1953). 

^ V. L. Ermolaev, E. N. Bodunov, E. B. Sveshnikova, and 
T. A. Shakhverdov, Nonradiative Energy Transfer of Elec- 
tronic Excitation (Nauka, Leningrad, 1977), in Russian. 

* V. M. Agranovich and M. D. Galanin, Electronic Ex- 
citation Energy Transfer in Condensed Matter (North- 
Holland, Amsterdam, 1982). 

^ Resonance Energy Transfer, edited by D. L. Andrews and 

A. A. Demidov (Wiley, Chichester, 1999). 
® V. L. Ermolaev, E. B. Sveshnikova, and E. N. Bodunov, 

Usp Fiz. Nauk 166, 279 (1996) [Phys. Uspekhi 39, 261 

(1996)]. 
■^ S. Jang, Y. J. Jung, R. J. Silbey, Chem. Phys. 275, 319 

(2002). 
® C. R. Kagan, C. B. Murray, and M. G. Bawendy, Phys. 

Rev. B 54, 8633 (1996). 
^ S. A. Crooker, J. A. HoUingd worth, S. Tretiak, and V. L. 

Klimov, Phys. Rev. Lett. 89, 186802 (2002). 
^^ R. Wargnier, A. V. Baranov, V. G. Maslov, V. Stsiapura, 

M. Artemyev, M. Pluot, A. Sukhanova, and I. Nabiev, 

Nano Lett. 4, 451 (2004). 
^^ A. Javier, C. Steven Yun, J. Sorena, and G. F. Strose, J. 

Phys. Chem. B 107 (2003). 
^^ A. Imamoglu, D. D. Awschalom, G. Burkard, D. P. Di- 

Vincenzo, D. Loss, M. Sherwin, and A. Small, Phys. Rev. 

Lett, 83, 4204 (1999). 
^^ E. Biolatti, R. C. lotti, P. Zanardi, and F. Rossi, Phys. 

Rev. Lett. 85, 5647 (2000). 
^* T. Holstein, S. K. Lyo, and R. Orbach, in Laser Spec- 
troscopy of Solids, edited by W. M. Yen and P. M. Selzer 

(Springer, Berlin - Heidelberg - New York, 1981), p. 39. 
^^ V. A. Malyshev, in Spectroscopy of Crystalls, edited by 

A. A. Kaplyanskii (Nauka, Leningrad, 1985), p. 100, in 

Russian. 
^^ T. T. Basiev, V. A. Malyshev, and A. K. Przhevusku, in 

Spectroscopy of Solids Containing Rare-Earth Ions, edited 

by A. A. Kaplyanskii and R. M. Macfarlane (North- 
Holland, Amsterdam, 1987), p. 275. 
^■^ S. Xia and P. A. Tanner, Phys. Rev. B 66, 214305 (2002). 
1* H. Sumi, J. Phys. Chem. B 103, 252 (1999). 
1^ K. Mukai, S. Abe, and H. Sumi, J. Phys. Chem. B 103, 

6096 (1999). 
^° G. D. Scholes, X. J. Jordanidis, and G. R. Fleming, J. 

Phys. Chem. B 105, 1640 (2001). 
^^ X. J. Jordanidis, G. D. Scholes, and G. R. Flemming, J. 

Phys. Chem. B 105, 1640 (2001). 
22 G. D. Scholes, Annu. Rev. Phys. Chem. 54, 57 (2003). 



^■^ S. Jang, M. D. Newton, R. J. Silbey, Phys. Rev. Lett. 92, 

218301 (2004). 
^-^ G. R. Fleming, G. D. Scholes, Nature 431, 256 (2004). 
'^■' A. P. Shreve, J. K. Trautman, H. A. Frank, T.G. Owens, 

and A. C. Albrecht, Biochim. Biophys. Acta 1058, 280 

(1991). 
2"^ R. Jimenez, S. N. Dikshit, S. E. Bradforth, G. R. Fleming, 

J. Phys. Chem. 100, 6825 (1996). 
2^ A. M. van Oijen, M. Ketelaars, J. Koler, T. J. Aartsma, 

and J. Schmidt, Science 285, 400 (1999). 
2® H.-M. Wu, M. Ratsep, R. Jankowiak, R. J. CogdeU, and 

G. J. Small, J. Phys. Chem. B 101, 7641 (1997). 
2® H. von Berlepsch, C. Bottcher, A. Quart, M. Regenbrecht, 

S. Akari, U. Keiderling, H. Schnablegger, S. Dahne, and S. 

Kirstein, Langmuir 16, 5908 (2000). 
*" H. von Berlepsch, C. Bottcher, A. Quart, C. Burger, S. 

Dahne, and S. Kirstein, J. Phys. Chem. B 104, 5255 

(2000). 
^^ H. von Berlepsch, S. Kirstein, R. Hania, C. Didraga, A. 

Puzglys, and C. Bottcher, J. Phys. Chem. B 107, 14176 

(2003). 
^2 C. Didraga, A. Pugzlys, P. R. Hania, H. von Berlepsch, K. 

Duppen, and J. Knoester, J. Phys. Chem. B 108, 14976 

(2004). 
^^ A. Pugzlys, P. R. Hania, C. Didraga, V. A. Malyshev, J. 

Knoester, and K. Duppen, in Ultafast Phenonena XIV, 

edited by T. Kobayashi, T. Okada, T. Kobayashi, K. A. 

Nelson, and S. De Silvestri (Springer, Berlin - Heidelberg 

- New York, 2005), Springer Series in Chemical Physics, 

vol. 79, p. 879. 
^"^ N. Fukutake, S. Takasaka, and T. Kobayashi, Chem. Phys. 

Lett. 361, 42 (2002). 
'^^ D. J. Heijs, V. A. Malyshev, and J. Knoester, Phys. Rev. 

Lett. 95, 177402 (2005); J. Chem. Phys. 123, 144507 

(2005). 
^•^ M. Bednarz, V.A. Malyshev, J. Knoester, J. Chem. Phys. 

117, 6200 (2002); 120, 3827 (2004). 
^"^ L. D. Landau and Lifshits, Quantum Mechanics (Perga- 

mon, Oxford, 1965). 
^® U. Weiss, Quantum Dissipative Systems (World Scientific, 

Singapore, 1993). 
^^ O. Kiihn and V. Sundstrom, J. Chem. Phys. B 107, 4154 

(1997). 
*° V. May and O. Kiihn, Charge and Energy Transfer Dy- 
namics in Molecular Systems (Wiley- VCH, Berlin, 2000). 
"^ T. Renger, V. May, and O. Kiihn, Phys. Rep. 343, 137 

(2001). 
^2 B. Briiggemann, K. Szenee, V. Novoderzhkin, R. van Gron- 

delle, and V. May, J. Phys. Chem. B 108, 13536 (2004). 
"^ M. Bednarz, V.A. Malyshev, J. Knoester, Phys. Rev. Lett. 

91, 217401 (2003). 



